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Abstract 

The initial value problem Vq, in all of the space, for the spatio 
- temporal FitzHugh - Nagumo equations is analyzed. When the re- 
action kinetics of the model can be outlined by means of piecewise 
linear approximations, then the solution of Vq is explicitly obtained. 
For periodic initial data are possible damped travelling waves and their 
speed of propagation is evaluated. The results imply applications also 
to the non linear case. 
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1 Introduction 



One of the reaction diffusion systems which models various important 
biological phenomena is given by 

+/(«) 

(1.1) 

[ || = bu -0v 

where the appropriate class of functions f(u) depends on the reaction kinet- 
ics of the model [1] - [4]. In the theory of nerve membranes, for example, the 
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system (1.1) is related to the spatio-temporal FitzHugh - Nagumo equations 
(FHN) with 



(1.2) 



f(u) — u(a — u) (u — 1) 



(0 < a < 1) 



and where u ( x, t ) models the transmembrane voltage of a nerve axon at 
distance x and time t, v(x, t) is an auxiliary variable that acts as recovery 
variable. Further, the diffusion coefficient e and the parameters b, (3 are all 
non negative [5], [9]. 

In addition to the propagation of nerve action potentials, (1.1) - (1.2) 
govern other several biological and biochemical phenomena; the list of ref- 
erences is very long and the variety of analytical aspects examined is wide. 
[10] - [15]. Further, we must remark travelling pulses and periodic wavetrains 
obtained by means of piecewise linear approximations of f(u) as 



where rj denotes the unit- step function. [16]- [21]. 

Typical boundary value problems related to the linear case (1.3) can be 
explicitly solved. Aim of this paper is the analysis of the initial value prob- 
lem Vo in all of the space; the fundamental solution and the explicit solution 
of Vq are determined when (1.3) holds. More, in the non linear case of the 
FHN model given by (1.1)-(1.2), the problem Vq is reduced to appropriate 
integral equations whose kernels are functions characterized by basic prop- 
erties. All this implies existence and uniqueness properties, together with a 
priori estimates. 



2 Statement of the problem and results 

Both the non linear source (1.2) and the linear approximation (1.3) in- 
volve a linear term — ku with k = a for (1.2), and k — 1 for (1.3). As 
consequence, the system (1.1) becomes 



(1.3) 



f(u) = r)(u-a) 



— u 



(0 < a < 1) 
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(2.1) 



u t — e u xx + ku+ v = f(u) 
v t + f3v - bu = , 



(x,t) e ft 



where f (u) = u 2 ( a + 1 — u ) when / is given by (1.2), while, in the 
linear case, ip (u) is equal to the constant fj that holds zero or one. 

The initial- value problem Vo related to (2.1) with <p (u) = fj is analyzed 
in the set 



(2.2) 



with the conditions 



tt T = {(x,t) : x eft, 0<t<T}, 



(2.3) 



u(x, 0) = uq(x), v (x, 0) = Vq(x), X e ft. 



When the linear approximation of / holds, then the explicit solution of 
the problem Vq can be obtained by means of functional transforms (Fourier 
with respect to x and Laplace as for t). If one puts formally 



(2.4) u(x,s) = / e~ st u(x,t) dt , v(x,s) = / e~ st v(x, t) dt , 

JO JO 



from (2.1) (2.3) one deduces: 



(2.5) u= ( [«o(0 + -] K [(x -£,s) de- / v (0Ki(x-Z,s) d£ 

JSR S J5R 



(2.6$ = b f MO + ^Kiix - £, s)di - 6 f v {i)K 2 {x - + ^4 

JK S JSR S + p 



where 
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with a 2 = s + k + -£3. 

s+p 

Now, let us consider the fundamental solution 



(2.8) 



9(x,y) 



2^ 

of the heat equation £ ^ — g y — k g = and, for n = 0, 1,2, let 



(2.9) G n (x,t) = j\-^g(x,y)(-^-) n7il J n _, (2y/b(t - y)y)dy, 

where J n {z) denotes the Bessel function of first kind. Then, if one puts 

(2.10) K = g(x,t) + G (x,t), K, = Gfat) (i = 1, 2), 

the following theorems hold. 

Theorem 2.1- In the half-plane 3?es > max( — k, —(3), the Laplace 
integrals of K n (x,t) (n = 0,1,2) converge absolutely for all \x\ > 0, and 
one has £ t K n (x,t) = K n (x,s). 

Theorem 2.2- The functions K , Ki, K 2 are C°° (Or) solutions of the 
integro differential equation : 

(2.11) z t -sz xx + kz + b / e _/3( *- r) z(x, r) dr = 

Jo 

and have the same basic properties of the fundamental solution (2.8) of the 
heat operator. Further, for i = andi =1, it results: 
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(2.12) Ki(x,t) = (d t +p)K i+1 , = (i = 0,1) 

while lim i^o * ) = only for \x\ > 0. 

These properties assure the convergence of the convolutions 

(2.13) K n *ij = j^K n (x-C,t)iP(OdC (n = 0,l,2) 
for all the functions that satisfy a growth condition of the form 



(2.14) 



[x)\ < c\ exp [c2 \x 



a+l 



, < a < 1 



with c\ and c 2 positive constants. 

Further, let consider the following functions 



(2.15)iV 1 (i) 



1 - e 2 COS (/y£) 



77(26 + /5 - /5 2 ) _iM t . 

■e 2 l sm(7t) 



2 7 (6 + /3) 



(2.16) N 2 (t) = 



1 — e 2 cos (7*) 



77(1 + ^) __i±£ t „ ; 



2 7 (6 + /3) 



e 2 sin( 7 £), 



where 7 = 
statements, t 



b - {13 - 1) 2 /4] 1/2 . Then, by (2.5)-(2.6) and the foregoing 
re explicit solution of the linear problem Vq is given by 



(2.17) 



u = Ni (t) + Mo * Kq — f * K\ 

< 

^ = N 2 (t) + M * Ki - Vq * K 2 + Vq (x) t~ fit , 

and the following conclusion is deduced. 
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Theorem 2.3 - When the data (uq, vq) are continuous functions that 
satisfy the growth condition (2.14), then the formulae (2.17) represent the 
unique solution of the problem Vq in the class of solutions compatible with 
(2.14). 



3 Travelling waves and a priori estimates 



A first example of applications is related to the linear case and concerns 
the analysis of travelling waves. By the explicit solution (2.17), for instance, 
when fj = 0, v , = u = Acos{w x), then one obtains: 

(3.1) u = (d t + /3)w, v = bw 

where 



, , A l+p+ew 2 . r . . 

(3.2) w = — e 2 sm{wx + at) — siii(to — at) 



2a 



with a 



2 



! . So, when b > ( 1+£ ^ 2 ~ /3 ) 2 , there exist 



damped travelling waves with speed equal to a/w. 

Moreover, in the non linear case, by (2. 5), (2. 6) one deduces the following 
integral equations 



u = u * K — v * K 1 + Jq (p(u) * K dr 

(3-3)4 

v — b uq * Ki — bvo * if 2 + v (x) e /3t + b f l ip(u) * K\ dr 

that imply a priori estimates. In fact, in the class of bounded solutions, if 
one puts 
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(3.4) ||m || = sup |mo(x)|, ||i>o|] = sup |vo(a;)|, ||<^|| = sup \(p(u)\ 

xG5R a;G3t «£SR 

by means of the basic properties of the kernels K , K\ K 2 , one obtains: 



(3.5) 



where 



\u\ < ci \\<p\\ + c 2 ( \\u \\t + \\v \\ ) E(t) 
\v\ < ci \ \(f\\ + c 2 ( ||u || + IKII * ) E(t) 



(3.6) 



E(t) = 



e -kt _ e -pt 
B-k 



and where the constants c±, c 2 depend on the parameters b, k, (3. Worthy of 
remark is the fact that the estimates (3.5) hold for all t, also when T — > oo. 
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